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Quantum information processing systems are often operated through time dependent controls; 
choosing these controls in a way that makes the resulting operation insensitive to variations in 
unknown or uncontrollable system parameters is an important prerequisite for obtaining high-fidelity 
gate operations. In this article we present a numerical method for constructing such robust control 
sequences for a quite general class of quantum information processing systems. As an application 
of the method we have designed a robust implementation of a phase-shift operation central to rare 
earth quantum computing, an ensemble quantum computing system proposed by Ohlsson et. al. 
P|. In this case the method has been used to obtain a high degree of insensitivity with respect to 
differences between ensemble members, but it is equally well suited for quantum computing with a 
single physical system. 
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Introduction 

Many potential quantum information processing sys- 
tems are controlled by means of a set of time-dependent 
parameters, such as quasi-static electromagnetic fields 
13 , radio-frequency, [3, , or optical fields 1 , 5 , 6 , . 
For most such systems, it is relatively simple to device 
a set of controls that implement a given evolution in an 
ideal situation. Often, however, a more careful choice of 
controls can lead to an implementation that is less sen- 
sitive to variations in unknown or uncontrollable system 
parameters. Examples of such robust implementations 
include system specific solutions such as the hot gate for 
ion trap quantum cornputing which is insensitive to vi- 
brational excitations [g, as well as more general tech- 
niques such as composite pulses, a technique originating 
in NMR spectroscopy 

In this article we describe a numerical method for de- 
signing robust controls for systems where the evolution is 
adequately described by a possibly non- unitary evolution 
operator U(t). This form does not allow a general mas- 
ter equation formulation, but it is sufficient to establish 
worst case behavior in many quantum computing settings 
where the worst case effects of decoherence and loss can 
be adequately modeled by a Schrodinger equation with a 
non-Hermitian Hamiltonian. 

As an application of the method, we will consider the 
construction of a robust phase shift operation for the rare 
earth quantum computing (REQC) system which 
is based on rare-earth ions embedded in a cryogenic crys- 
tal. In each ion, two metastable ground-state hyperfine 
levels, labeled |0) and serve as a qubit register which 
is manipulated via optical transitions from both states to 
an inhomogeneously shifted excited state |e). The REQC 
system is an ensemble quantum computing system and 
macroscopic numbers of ions are manipulated in paral- 
lel, addressed by the value of the inhomogeneous shift 
of their |e)-state. To obtain a sufficient number of ions 
within each frequency channel, it is necessary to operate 



on all ions within a finite range of inhomogeneous shifts, 
and the main difficulty in operating the REQC system 
is to achieve the same evolution for each of these ions 
independent of their particular inhomogeneous shift. 

This article is divided into two sections: in section^we 
describe the method we have used to design robust gate 
implementations; these results should be applicable to a 
variety of quantum information processing systems. In 
section ^ we present the results of applying the method 
to a specific problem relating to the REQC system, and 
show that it is indeed possible to obtain very high degrees 
of robustness. 



I. DESIGNING ROBUST GATE 
IMPLEMENTATIONS 

We consider a collection of quantum system which 
evolve according to a set of time-dependent controls 
£{t). In addition to e, the single system Hamiltonian 
-ff (^, £(i)) depends on a system specific set of uncontrol- 
lable or unknown parameters ^, such as field strength 
or quantum numbers corresponding to unused degrees of 
freedom. 

The evolution of each system is governed by the 
Schrodinger equation, 

zft^U(t) = H{t e{t)) U(t), U(0) = 1, (1) 

where we will allow the Hamiltonian H to include non- 
Hermitian terms describing loss and decoherence. 

Our goal is to choose a set of controls that lead to an 
evolution U(r) which is as close as possible to a given de- 
sired evolution Uq over a range of ^-values. To quantify 
this, we introduce an objective functional J(^, e) which 
describes the performance of a set of controls e for a given 
value of ^. By convention we take a low value of J to in- 
dicate a good performance, and the problem of finding a 
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robust set of controls thus corresponds to minimizing 

Jx(e) =niaxJ(^,£), (2) 

where X is the set of ^-values for which we want the im- 
plementation to perform well. The conceptually simple 
approach we have taken to this problem is to replace X 
with a discrete subset X' <Z X , so that the minimiza- 
tion of Jx' has the form of a standard minimax problem, 
which can be solved efficiently provided that we are able 
to calculate dJx' /de. Below we show how to achieve this 
by methods from optimal control theory. 



A. Calculating dJ/de 

In this section we show how dJjde may be calculated 
for a quite general class of objective functionals. To 
keep the notation simple and avoid unnecessary restric- 
tions, we will consider the following generalization of the 
Schrodinger equation (^), 



U(i) -/(U(t),£(t)), U(0)-1, 



(3) 



determining the evolution of a complex- valued, time- 
dependent matrix U due to a set of real-valued, time- 
dependent controls, sit). We will consider objective func- 
tionals of the form 



J(e) = 0(U(T)) 



?(U(t),e(i))dt, 



(4) 



where 0(U(T)) is a real-valued function quantifying how 
close the final state U(T) is to our goal, and the real 
valued function Z((U(i), e(t)), referred to as a penalty 
function, can be chosen to discourage the use of certain 
control values. Our goal is to calculate dJ/de subject to 
the constraint that e and U obey Eq. To achieve 

this we introduce the modified objective functional: 



Trf At 



U-/(U,£) +\i.c.]dt, (5) 



J' = J- 



where the complex time-dependent adjoint state matrix 
A is in effect a continuous set of Lagrange multipliers 
leaving J' identical to J, provided that U and e obey 
Eq. If we require A to obey the adjoint equations 



At 



dh 



At(T) = 



(6) 



where U and IJt should be considered as independent 
with respect to the partial derivative and h is defined as 



h{e, U, A) = Z(U, £) + Tr (At/(U, e) + h.c. 



(7) 



we find by integration by parts that the differential of J' 
is given by [ll| 

t in 

(8) 



dJ' ^ I ^ 5e{t) dt, 
Jo de 



from which the derivatives of J with respect to the pa- 
rameters used to parametrize £ can be calculated by the 
chain rule. 

We now return to the case of a quantum system gov- 
erned by the Schrodinger equation If we assume the 
penalty function I to be independent of U, the adjoint 
equations in this case are 



d 



^n-m - Wit sit)) Ait), A(T) = ^ 
and dJ is given by JHJl with 



dh dl 2 f.idH^ 

ir = 7r + -;^iinTr At— u 



(9) 



(10) 



The role of the adjoint state and the adjoint equations is 
often described as back-propagating the errors in achiev- 
ing the desired final state. If H is Hermitian, the bound- 
ary value for A can be optimized for numerical compu- 
tation as shown in Appendix IbI 



B. Fidelity of quantum evolution 

We will now discuss the choice of the function 0, quan- 
tifying how well the obtained evolution U(T) approxi- 
mates Up. As we are concerned with quantum informa- 
tion processing we will assume that all operations start 
out with an unknown state in the qubit subspace Tiq of 
the full system Hilbert space Ti, and that this subspace 
is invariant under the ideal evolution Uq. The function 
should not depend on the evolution of states outside Ti-q, 
nor on collective phases on the states originating in Uq. 
A cautious choice of (/) fulfilling these conditions could be 
based on the worst case overlap fidelity jl2|: 



T — min 



utu(T) 



(11) 



which measures the least possible overlap between the ob- 
tained output state U(T) j?/'} and the ideal output Uq IV') 
for initial states in TLq . This fidelity measure has the de- 
sirable quality that both population transfer from Tiq to 
TLq and population transfer completely out of 7i, as de- 
scribed by a non-unitary evolution, is counted as loss of 
fidelity. 

From the point of view of optimal control, a significant 
drawback of the worst case overlap fidelity, T , is that it 
is computationally complicated |l3| . A computationally 
accessible fidelity measure which share many appealing 
features with T is the trace fidelity 0], 



(12) 



where n = dim (Tig). As shown in Appendix 1X1 T is 
related to J- by the strict bound 











Tr (\ 


jtu) 


n 







l-J^<n(l-r), 



(13) 
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indicating that we can safely replace .7-" by T for numer- 
ical computations on a few qubits at high fidelity. 

For numerical calculations it is beneficial to use 
rather than T in the calculations presented in the 
next section we have used (f> = \ — T^, leading to an 
adjoint state boundary condition of 



A(T) 



d(j) 
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1 



-^Uo Tr(USU(r)), 



(14) 



which can be directly computed. 



C. Minimization algorithms 

One approach to minimizing J is to directly solve the 
extremum condition dJ /de — for e. This task is sig- 
nificantly simplified if a penalty function proportional to 
£^ is introduced, h = + ho, so that the extremum 
condition according to ((HJ reads, 



eit) 



1 dhp 



(15) 



which may be used as an iterative formula for calculat- 
ing e. Variations over this iterative app roach give rise to 
the Krotov and Zhu-Rabitz p^ll8l | algorithms which 
have been shown to have excellent convergence proper- 
ties, and have been successfully applied to optimal con- 
trol of unitary transformations for one set of parameters 
by the group of R. Kosloff [H H^. A unifying view of 
these direct methods can be found in Ref. [13 . 

In the work presented here we have chosen to use 
an indirect minimization algorithm: rather than trying 
to solve the extremum condition directly, we have used 
the gradient information obtained through Eq. JSJ as 
input for a general sequential quadratic programming 
procedure based on a constrained quasi-Newton method 
[20I I21I |23 | . The primary advantage of this approach is 
that we have total freedom to choose the parametriza- 
tion of the controls, and can place arbitrary bounds on 
these. This allows us to more accurately model the fact 
that the experimental limitations most often only distin- 
guish between possible and impossible controls: no pos- 
sible controls are significantly harder than others. An 
explicit field strength limit also serves to introduce an 
absolute scale on which to introduce decay strengths etc. 



II. APPLICATION TO RARE EARTH 
QUANTUM COMPUTING 

The motivation for the work presented in this article 
has been the design of robust gate implementations for 
the REQC system mentioned in the introduction. As an 
example, we will consider the construction of a robust 
implementation of the single qubit operation 



which could simply be implemented by a single 27r pulse 
on the |0)-|e) transition if we were not concerned with 
robustness. Our primary concern will be to make the 
implementation robust with respect to variations in the 
inhomogeneous shift S of the |e)-state in order to allow 
the use of finite with channels. Since it is experimentally 
difficult to obtain a homogeneous field strength over the 
crystal, we would also prefer the implementation to be in- 
sensitive to variations in the relative field strength, which 
we will denote 7. 

In addition to requiring the implementation of Uq to be 
robust with respect to variations in S and 7, we will add 
the requirement that ions outside the channels should 
not be affected, as this allows us to use the obtained 
implementation of Uq as a part of a controlled phase 
shift operation 0,0: If the |e)-states of the qubit ion 
and a controlling ion are coupled sufficiently strongly by 
static dipole interaction, an excitation of the controlling 
ion will effectively shift the qubit ion out of the channel 
thus conditioning the evolution of the qubit on the state 
of the controlling ion. 

Even though the simplest implementation of Uq for 
an ideal ion with S = and 7 = 1 would involve only 
the |0} and |e) levels, a robust implementation must also 
involve the |l)-state since the coupling of |0) to |e) will 
result in a 5-dependent phase on |0} which can only be 
compensated by introducing the same phase on the |1) - 
state, e.g. through phase compensating rotations [l3l| . 
A highly successful example of this approach is the sech- 
pulse sequence suggested by Roos and M0lmer J3] , which 
as illustrated by Fig. ^a) yields a very robust imple- 
mentation, achieving high fidelities over a wide range of 
parameter values. 

We model the REQC system by the single ion Hamil- 
tonian 



H 



-S\e) (e| 



^ (17,(0 |e)(z|+h.c.), (17) 

1=0,1 



Uo = |l)(l| 



(16) 



which does not include any effects of decay or decoher- 
ence. In the notation introduced in section ^ the system 
parameters are ^ = (7, S), the controls are e = (JIq, i^i), 
and the qubit subspace is Hq = {|0) , |1)}. No penalty 
function is used: we use J = 4> = 1 — T^, and limit the 
field by strict bounds on fli{t), as this is the relevant 
limiting parameter in the REQC system. Inspired by the 
success of the sech-based solution and the hat-like Fourier 
spectrum of the sech-pulse, ili{t) is parametrized in terms 
of a truncated Fourier basis. Based on trial and error we 
have arrived at 49 ^-values to constitute X', some within 
the neighborhood of (5 = 0, 7 = 1, and some at large 
detunings where the ions should not be disturbed. The 
result of the optimization with this choice of X' is shown 
in Fig. n^b), where the circles indicate the members of 
A'. It is evident from the plot that the optimization has 
achieved a high fidelity over an even larger range of pa- 
rameters than the sech-pulse sequence illustrated in Fig. 
da). 

With respect to not disturbing the detuned ions, both 
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(a) sech-pulses 



(b) Optimized pulse 



FIG. 1: The worst case overlap fidelity of Uo as given by Eq. I16II implemented (a) by a series of sech pulses as suggested 
by Roos and M0lmer j2^ . and (b) by an optimized pulse. The fidelity is plotted as a function of relative field strength 7 and 
inhomogeneous shift 5 relative to the maximal resonant Rabi frequency ilo- The duration of both pulses is 247r/r2o, and the 
circles in (b) indicates the parameter set X' used during the optimization. 




FIG. 2: The effect of the sech-pulse sequence (dashed line) 
and the optimized (solid line) implementation of Uo on far- 
detuned ions, as described hy 1 — where is calculated 
with respect to the identity, as these ions should ideally not be 
disturbed. Both implementations achieve worst case fidelities 
very close to unity, with the sech-pulse sequence achieving the 
best results. The curve for the optimized pulse is actually a 
running maximum, as the actual value of 1 — oscillates with 
a period of 2n/T. 



III. CONCLUSIONS AND OUTLOOK 

We have shown that it is possible to construct highly 
robust gate implementations for quantum information 
processing by a quite general method. In particular, 
the method has been used to greatly enhance the per- 
formance of a gate implementation for a model REQC 
system by extending the range of inhomogeneous shifts 
and relative field strengths over which an acceptable per- 
formance is achieved. 

The model REQC system used in this article ignores 
many performance degrading factors, the two most im- 
portant being decoherence and implementation noise. 
Decoherence could in the present case be adequately 
modeled by a non-Hermitian Hamiltonian, for which we 
expect the method described in this paper to be able to 
find a robust implementation as in the decoherence-free 
case. It is not clear, however, how the method could 
be extended to address the problem of robustness with 
respect to implementation imperfections. 



the optimized pulse and the sech-pulses obtain fidelities 
within 10^^ of unity for \S\ > 5Q,o, where Qq is the max- 
imal resonant Rabi frequency at 7 = 1. As illustrated by 
Fig. 13 which only shows fidelities at 7 = 1 as U(T) is 
nearly independent of 7 at \S\ ^ fio, the sech-pulse se- 
quence performs better than the optimized pulse in this 
regime. 
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APPENDIX A: THE TRACE FIDELITY 

In this section we prove the relation 1)13(1 between the 
worst case overlap fidelity, T , and the trace fidelity T. 
Referring to the definition Hll|l . we note that T is com- 
pletely determined by the restriction O of the operator 
UjU(r) to Uq. Since U(T) describes the evolution of 
a quantum system, it is possible to extend it to a uni- 
tary operation on a Hilbert space containing 7Y, and O is 
consequently the restriction of a unitary operator to Tig . 
In the ideal case O will be equal to the identity on Uq, 
perhaps with the exception of a complex phase. 

T is defined as the minimum of the overlap | ("0101?/;) | . 
Since the unit sphere of C" is compact, this minimum will 
be attained for some T = \ (V'o|C)|'0o) I- We now 

extend HV'o)} to an orthonormal basis {|V''fc)}fc=o....,n-i 
by the Gram-Schmidt process. Evaluating the trace fi- 
delity in this basis we find by the triangle inequality: 



n — 1 

r < -y^\{i>H\o\i,K)\ 

n ^ — ^ 



(Ala) 
(Alb) 



where we have used that | (-i/jlOl-!/)) | < 1 for all since 
O is the restriction of a unitary operator. By rewriting 
(|Alb|l we obtain the desired relation, 1 — < n(l — T). 

We note that the established bound is strict in the 
sense that for any < J-q < 1, the operator 

0^ = l-(l-J-o)|V)(^| (A2) 
will fulfill Eq. (|13|l with equality for any \tp). 



APPENDIX B: OPTIMIZED ADJOINT STATE 
BOUNDARY CONDITION 

In the case of a Hermitian Hamiltonian and a penalty 
function I, that does not depend on the state U, it is 
possible to modify the adjoint state boundary condition 
(fnji to reduce the required accuracy of the adjoint state 
propagation. 



In this case, we find according to Eqs. © and H10() 
that dJ has the form 



dJ 



2Im 



dH 



Se{t)dt, (Bl) 



where Xk and Xk denotes the fc-th columns of U and A 
respectively. Since H and thus dH/de are assumed to be 
Hermitian, dJ as given by HB1|) is not affected by adding 
to Xk a component of akXk for any real ak- Furthermore, 
since A^. and Xk evolve according to the same Schrodinger 
equation, this corresponds to modifying the boundary 
condition for the adjoint state to read: 

XkiT) = -^+akXkiT), (B2) 

for any real ak- 

The obvious use of the freedom in the choice of bound- 
ary value is to minimize the norm of the adjoint state, in 
order to relax the requirements of the relative precision 
of the adjoint state propagation. This minimum is easily 
calculated from (|B2|) , but we prefer to illustrate the phys- 
ical background of the result by calculating it in a dif- 
ferent way: The freedom in the choice of boundary value 
(|B2|) is allowed by the Hermiticity of H. The same Her- 
miticity ensures that Xk is normalized, so that (t)on{xk), 
where n is the normalization function n{v) = v/\v\, is 
equal to 4>{xk). The gradient of (/) = o n, however, is 
different from that of 6. In fact we find 



dxl 



dxl 



Re 



'dxl 



Xk-, 



(B3) 



where Xk and x\, should be considered independent with 
respect to the derivative. Comparing this expression to 

Eq. I)B2)I . it is tempting to let Uk = ^e{x\d(f> / dx\) , 

which is indeed the answer found by minimizing |a;fc(T)| 
subject to Eq. (|B2)| . 

The modified boundary condition of Eq. ((B3)l carries 
the same error information as that of Eq. (|14|l . but the 
required relative numerical precision when propagating 
the adjoint state will be significantly reduced. 
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